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Abstract 
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pendent or independent of brain size. It is discussed how this formula 
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1 Introduction 



The salient feature of a macroscopic cortical organization is the presence of different 
functional units such as columns and areas. These units form an ordered network 
of connections mediated by axonal bundles in white matter. On a microscopic scale, 
however, neurons are connected in sparse, stochastic local circuits (Braitenberg, 1978a; 
Douglas and Martin, 1991). The transition between the two levels of organization takes 
place on a column-size scale. 

It has been argued that both microscopic and macroscopic connectivity, i.e., the 
fraction of connected sites, should decrease with brain size (Stevens, 1989; Ringo, 1991; 
Ringo et al, 1994), because this would decrease the total axonal length; a feature de- 
sirable especially in bigger brains. This expectation follows from the so-called min- 
imal axon length principle, which assumes that the brain tries to save its biochemi- 
cal resources (Mitchison, 1992; Cajal, 1995; Cherniak, 1995; Murre and Sturdy, 1995; 
Chklovskii and Stevens, 1999). 

The purpose of this article is to investigate whether this is a justified presumption by 
studying how both microscopic and macroscopic connectivity depend on brain size. It is 
found that, although the connectivity between neurons decays with brain size (a known 
fact), the connectivity between cortical areas is either weakly decreasing or invariant 
with the brain volume. We argue that the formula for the connectivity between areas 
can be used to estimate an average axon length in white matter, thus providing a useful 
practical tool, especially in the face of lacking data. We also suggest that brains, in order 
to perform efficient computation, have to sacrifice some of their computational power, 
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and this is due to a trade-off between managing brain size and its limited biochemical 
resources, and maintaining functional operativeness. We speculate that this trade-off 
may lead to observed scaling laws between some cortical parameters. 

In all considerations below, it is assumed that the brain volume Vf, and the gray 
matter volume V g scale with an exponent close to 1, which is well justified experimen- 
tally (Jerison, 1973; Frahm et al, 1982; Prothero and Sundsten, 1984; Hofman, 1985; 
Hofman, 1989; Barton and Harvey, 2000). Therefore, we use terms "scale with gray 
matter volume" and "scale with brain volume/size" interchangeably and treat them as 
equivalent. 

For completeness, first let us determine how neuronal connectivity depends on brain 
size. An average connectivity p (or probability of connection) between neurons can be 
defined as p = M/N, where M is the average number of synapses per neuron, and N is 
the total number of neurons in the cortex. Since the volume density of synapses in gray 
matter is brain size independent (Schiiz and Demianenko, 1995; Braitenberg and Schiiz, 
1998), i.e. NM/Vg = const, then we obtain that M ~ V g /N, and as a consequence 
p ~ Vg/N 2 . The total number of neurons N is proportional to the total cortical surface 
area W (Rockel et al., 1980). The latter scales with the brain volume as: W ~ V g ' 9 
for large convoluted brains (Hofman, 1985). This leads to the following scaling between 
the average connectivity p and brain size for convoluted brains 

P-v- - 8 . (1.1) 



3 



Thus p decreases quickly with brain size. As an example, the human brain volume and 
the rat brain volume differ by a factor of 614.6 (Hofman, 1985). From this, it follows 
that the average connectivity p is about 170 times smaller in human than in rat. 

In the next Sections we study scaling of the area connectedness with the brain 
volume. 

List of symbols used in this article 

a is a dimensionless parameter characterizing cortical geometry and a pattern of axonal 
organization in white matter, 

c is the ratio of the volumes of gray matter and the whole brain, 

d is the average diameter of axons in white matter, 

E is the basal metabolic energy rate used by gray matter at rest, 

/ is the fraction of active excitatory synapses in gray matter at resting conditions, 

K is the number of areas in the cortex, 

L is the average length of axons in white matter, 

M is the average number of synapses per neuron in the cortex, 

iV is the total number of neurons in the cortex, 

< n > is the average degree of separation between cortical areas, 

p is the average probability of connection between neurons in the cortex, 

q is the probability of sending at least one macroscopic axonal bundle to white matter 

by a module, 

Q is the average probability of connection between two areas, 

P n is the probability that two areas are connected in at least n steps, 
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S n is the probability of connection of two areas via at least one of the paths that uses 

n intermediate areas (steps), 

W is the total cortical surface area, 

W is the surface area of one cortical area, 

V g is the gray matter volume, 

V w is the white matter volume, 

Vb is the brain volume, 

Vbody is the volume of the whole body, 

a is the scaling exponent of the cortical area number with the gray matter volume, 
(5 is the scaling exponent of the axonal length in white matter with the gray matter 
volume, 

7 is the scaling exponent of the white matter vs gray matter volume, 

5 is the scaling exponent of the cortical area connectedness with the gray matter volume, 

k, is the probability of connection between a given module in one area to some other 

area, 

k is a fraction of connected cortex by a module, 

£ is the linear size of a module in the cortex, 

r is the average conduction delay between cortical areas. 

2 Connectivity between cortical areas 

In this Section a probability of a connection between cortical areas is derived. We 
assume that the cerebral cortex of the total surface area W is divided into K areas of 
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an identical surface area W . Each of the areas, in turn, is composed of modules of 
linear size £. We define a module as a local group of neurons with similar functional 
properties that is capable of sending at least one macroscopic (containing at least several 
axons) coherent bundle of long-range (cortico-cortical) axons to a particular place in 
the cortex (Pandya and Yeterian, 1985; Zeki and Shipp, 1988; Felleman and Van Essen, 
1991). One can think about such defined modules as being cortical columns or barrels 
that have been found in the visual and somato-sensory cortices of different mammals. 

Consider two cortical areas A and B. We assume that the area A is connected to the 
area B if at least one of the modules in A is connected to B. Note that this definition 
of connectivity ignores the issue of strength of the connection, and focuses only on 
the existence of a link between areas. Also, we neglect the explicit influence of the 
distance between areas on their connectivity, since no precise quantitative anatomical 
data exists. In this sense, we consider only average area connectedness. However, we 
include implicitly the spatial effect by considering the average length of axons connecting 
different areas. If we assume that there are limited biochemical resources in the cortex, 
then in particular, the length of axons that modules can send to white matter should 
be kept as small as possible. This restriction puts a spatial constraint on the number 
of areas a module can reach, and this in turn, influences the probability of connection 
between areas. An average probability k that a given module in A connects with any 
module in B (probability of connection between a module in A and the whole area B) 
is given by 
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K = q(W /W)(aL 2 /W), 



(2.1) 



where q is the probability that the module in A sends at least one macroscopic axonal 
bundle to white matter. This probability should be close to 1, since e.g. for the mouse 
cortex pyramidal cells, many of which project to white matter, constitute 85% of the 
total number of cortical cells (Braitenberg and Schiiz, 1998) and a module consists 
of about 10 3 neurons (Braitenberg, 1978b). Thus, it is very likely that a sufficient 
number of axons is sent to form at least one coherent bundle, and this process should 
not depend on brain size. The ratio W /W in eq. (2.1) is the probability that this 
bundle terminates in the area B. In general, this probability can be different for a 
different pair of cortical areas (c.f. Young et al, 1995) and does not have to be equal 
to Wo/W = 1/K. However, since we are interested only in average probability, the 
assumption of uniformity is reasonable and sufficient. The dimensionless ratio aL\jW 
is the probability that the area B can be physically reached given finite average axonal 
length Lq in white matter (fraction of areas that can be connected given finite fiber 
range). The factor a is some dimensionless constant characterizing a particular cortical 
geometry (either convoluted or smooth) and a pattern of axonal organization in white 
matter (e.g., how many, on average, bundles a module sends to white matter). This 
factor can be weakly species-dependent, however, it should not depend in any systematic 
manner on brain size. We do not make any specific assumptions about the cortical or 
axonal geometry, and therefore our results and conclusions drawn are quite general. 
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The latter do not depend on a particular value of the parameter a. 

The probability Q that the area A connects with the area B (area connectedness) 
is exactly complementary to the probability that none of the modules in A connects 
with the area B, i.e. (1 — k) w °^ 2 , where W /£, 2 is the number of modules in every area. 
Thus, our Q is given by 

Q = 1-(1-k) w ^ 2 . (2.2) 

We can simplify this formula by using an identity 1 — k — exp[ln(l — k)], and by 
expanding the logarithm in it for small k according to ln(l — k) ~ —k. This produces: 
(1 — k) w °^ 2 ps exp(— kWo/C, 2 ). The next step is to substitute the right-hand side of 
eq. (2.1) for k, and to use the fact that W = W/K. After these operations, we finally 
obtain the following expression for the average connectivity Q: 

e~l-<*p(-f§). (2.3) 

From this formula, which is one of the main ones in this paper, it follows that the 
average area connectedness for a given cortex depends on its four basic characteristics: 
cortical geometry, the average length of axons in white matter, module size, and the 
number of cortical areas. 

It is also possible to find an expected number of modules in one area that connect 
with another area. Assuming that modules are statistically independent, i.e. the prob- 
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ability of sending axonal bundles for a given module does not depend on other modules, 
the distribution of the number of modules in A reaching B is represented by a binomial 
distribution. Thus the average number of modules in A connecting with B is given by 
the product of the probability that a module in A connects with area B (k) and the 
number of modules in A (Wo/!; 2 ), i.e. aqL^/^K 2 . The latter expression is exactly 
the same as the argument in the exponent in eq. (2.3), and this fact provides a useful 
interpretation of that equation. 

To find how Q scales with brain size, one therefore first has to determine how the 
expected number of connected modules between two areas depends on the brain volume. 
This is equivalent to finding how the above basic characteristics scale with brain size. 
This will be done in the next section. Here, we estimate the average connectivity Q 
for the mouse cortex, and also, the average length of long-range axons for few other 
species. All the numbers provided are for one isolated hemisphere only. 

Mouse is the only animal for which the distribution of axonal length in white matter 
of one hemisphere has been measured directly (Greilich, 1984; Braitenberg and Schiiz, 
1998, Fig. 62). From this, one can estimate that L ~ 3 mm. Module size £ can be 
assumed to be of the order of an average barrel size, i.e. £ 2 ~ 0.1 mm 2 . The unknown 
factor a in eq. (2.3) can be estimated from data of Schiiz and Liewald (Schiiz and 
Liewald, 2001). They found, using injections of anterograde tracer BDA, that neurons 
under a surface area of about 0.1 mm 2 , i.e. of approximately a barrel size, in the 
mouse cortex project onto 12% of the total cortical surface area of one hemisphere. 
The corresponding number for the macaque cortex is about 1 — 2% (Schiiz, private 
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communication). To find a we use eq. (2.1) for the probability of connection between 
a module in A and the area B. In the present case, however, the surface area of B 
is actually the whole cortex (of one hemisphere), i.e. W = W , and we obtain the 
following equation: 

R = aqLl/W, (2.4) 

where R = 0.12 is the fraction of mouse cortex which is connected by a module. The 
total cortical surface area is a ratio of the gray matter volume, which in one hemisphere 
is V g = 56 mm 3 , and cortical thickness, which is about 0.85 mm (Braitenberg and Schiiz, 
1998). Combining all the numbers, we obtain a xs 0.87 (for which can be used 

to determine the area connectedness. The number of cortical areas for mouse depends 
on the methodology used (see below). If we use the Caviness criteria (Caviness, 1975) 
then K = 24. However, if we use methodology of Kaas et al, then K = 6 — 8 (Krubitzer 
and Huffman, 2000). In the former case, this leads to Q — 0.13, while in the latter this 
yields Q = 0.71 - 0.89. 

For a few other mammals, i.e. rat, cat, and macaque monkey, it has been possible 
to determine, based on experimental data, the average connectivity Q. Unfortunately, 
we are unable to compare those numbers with predictions of eq. (2.3) due to the lack 
of data for L in these species. However, we can perform a reversed computation and 
evaluate the axonal length for each of these animals directly form eq. (2.3) [i.e. solve 
eq. (2.3) for L ], that is 
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This formula, as well as formula (2.4) can have practical importance for estimation 
of average length of fibers in white matter, since it is technically easier to determine 
the connectivity Q, or the fraction k of connected cortex, than to measure directly L . 
However, the area and connectivity counts, must be done within a single methodology 
(see below). 

For all the above animals the factor a and the sizes of modules are roughly equal 
(this will be justified in the next section) to the corresponding numbers for mouse. 
Additionally, for rat we can find the number of connections between areas based on 
a connectivity graph in Kolb (Kolb, 1990; p. 26, Fig. 2.1). The total number of 
connections between areas (we assumed that all of them are reciprocal, which is a 
reasonable assumption, although there is no information about this in the article) is 
found to be 168, which gives 6.46 connections per area. Since the number of areas 
is K = 26 (Kolb, 1990), we estimate that the connectivity for the rat cortex in one 
hemisphere is Q = 6.46/26 = 0.25. Using eq. (2.5), this yields L = 4.7 mm. For 
the cat cortex with Q = 0.27 and K = 65 in one hemisphere (Scannell and Young, 
1993; Scannell et al, 1995; Young et al, 1995), we obtain L = 12 A mm. Finally, for 
macaque with Q = 0.15 and K = 73 (Young, 1992; Young, 1993; Young et al, 1995), 
we obtain axonal length L = 10.0 mm. The latter value, which is somewhat smaller 
than expected may result from, for instance, underestimation of the connectivity Q 



aq Kl-Q, 



(2.5) 
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(Young et al, 1995, p. 130) or the parameter a. However, we can perform a second, 
independent, evaluation of L using eq. (2.4). Taking the macaque cortical surface area 
of one hemisphere as W — 1.25 • 10 4 mm 2 (Hofman, 1985; Stephan et al, 1981), we find 
L = 12.0 — 17.0 mm (for k = 0.01 — 0.02; Schiiz, private communication), which seems 
to be more reasonable. 

3 Scaling of basic cortical characteristics and area 
connectedness with brain size 

In this Section we study how the number of cortical areas and their size, module size, 
the long-range axon length, and the area connectedness depend on the brain volume. 

First, let us consider the scaling of the number K of cortical areas with brain size. 
Although, there is no consensus on how to define areas, and different physiological 
criteria lead to different counts, it is possible to obtain such a scaling law if one focuses 
on a single methodology. This has been done by Changizi (Changizi, 2001), who used 
the methodology of Kaas and Krubitzer (Kaas, 1987; Krubitzer, 1995; Krubitzer et al, 
1997). He found for 11 species with brain volumes spanning 2 orders of magnitude 
that K ~ V™ with a ~ 0.4. It is rather unlikely that another methodology would 
change this exponent drastically. It is also interesting to note that a similar scaling 
law has been proposed theoretically for hypothetical cortical compartments forming a 
completely connected network (Braitenberg 1978b; Braitenberg, 2001). 

We can combine the above scaling law for the number of cortical areas K with the 
previously presented (in Introduction) scaling law for the total cortical surface area W 
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to find a useful allometric relation between the average cortical area size Wq and the 
brain volume. Since W ~ V g 09 (Hofman, 1985), and K ~ V g 0A (Changizi, 2001), and 
because W = W/K, we obtain W ~ V^ - 5 . Thus, the size of cortical areas increases 
moderately with brain size. 

The dependence of a module size on the brain volume can be determined by noticing 
that the module size should correlate strongly with the size of cortical patches. The 
dependence of the size of patches in VI of different species on the gray matter volume 
V g is shown in Fig. 1. This figure reveals that there is no systematic dependence of £ 
on V g , and this suggests that £ is invariant with respect to brain size. 

As was mentioned above, there are no data on the average length of long-range 
axons in mammals (with mouse being an exception), and therefore we have to make a 
reasonable assumption about its scaling with brain size. If we assume that the average 
axon length L in white matter is proportional to its diameter, then L ~ V^ 3 , where 
V w is the white matter volume. The white matter volume is related to the gray matter 
volume via an allometric scaling law: V w ~ VJ, with 7 = 1.22 — 1.33 (Frahm et al, 
1982; Hofman, 1989; Allman, 1999; Zhang and Sejnowski, 2000; Barton and Harvey, 
2000). It seems that the difference between 1.22 and 1.33 in the exponent 7 results 
from the fact that different researchers included different species in their analysis. For 
example, the white matter volume of insectivores scales against the rest of the brain 
volume with a lower exponent than a corresponding exponent for primates (Barton and 
Harvey, 2000). Combining the above, we obtain for the average axon length in white 
matter L ~ ~ VJ^ 3 = V@, where (3 = 7/3. Thus the exponent (3 is in the range 
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0.41 — 0.44. We can verify this result independently. In Fig. 2 we plot our values of 
Lq, estimated in Sec. 2 based on eqs. (2.4) and (2.5), against the gray matter volume. 
The figure yields the value of the exponent (3 about 0.26, which is clearly different from 
the value ~ 0.4 above. The source of the discrepancy between the two numbers may be 
twofold. First, the assumption about the linear relationship between axonal length in 
white matter and the white matter diameter may not be quite correct. Second, our Fig. 
2 uses only 4 mammals and that may be a reason for a poor statistics. It is possible 
that the true value of the exponent (3 lies in the range 0.26 — 0.44. 

Now we are in a position to determine how the connectivity between cortical areas 
scales with brain size. Combining the scaling laws for L , £, and K, we find that the 
expected number of connected modules between two areas aqLl/^K) 2 scales as 

a 1 L l = A T/2(/3-a) 

^K 2 9 ' 

from which we obtain the scaling relation for the area connectedness 

q „ i _ exp (-AV- 5 ) , (3.2) 

where 5 = 2(a — (3) and A is some positive brain size independent constant. It is 
apparent that for the scaling of Q with brain size, it matters only the relative scaling 
of the average axon length and the number of areas. In the case when (3 ~ 0.4 the 
value of the exponent 5 ~ 0, suggesting almost perfect independence of Q of brain size. 
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(3.1) 



However, in the second case when (3 = 0.26, there should be a slow decay of Q with the 
brain volume, since then 5 = 0.28. Below, we estimate the magnitude of this decay. 

The constant A appearing in eqs. (3.1) and (3.2) can be determined by using our 
previous estimates of a = 0.87 and £ 2 = 0.1 mm 2 , data in Fig. 2 on L scaling, and 
data on K scaling (Changizi, 2001). From Fig. 2 we obtain that L = 1.10 V®' 26 . From 
Fig. 3 in Changizi (Changizi, 2001) we find the area number K = 0.42 V b a4 , where 
14 is the brain volume of one hemisphere expressed in mm 3 (note that Changizi uses 
the total brain volume of two hemispheres). However, the cortical area count used by 
Changizi (Kaas-Krubitzer methodology) underestimates the area number by a factor 
of 2.6-4.0 (calculation for monkey, cat, and mouse) in comparison to more traditional 
approaches (Caviness, 1975; Kolb, 1990; Young et al, 1995). In order to be consistent 
with the previous sections, where we used those traditional approaches in our estimates 
of the axonal length, we rescale the Changizi scaling law by a factor of 3.3, which is 
the average of the above numbers. Thus, a modified expression for the K scaling takes 
the form: K m 3.3K w 1.45 V b ' 4 , where K is the modified cortical area number. We 
also need the relationship between the brain volume and its gray matter volume 
V g that takes the form: V g /Vt, = c, where c is approximately constant within a single 
mammalian group. For example, c is about 0.5 for primates and 0.1 for insectivores 
(Stephan et al, 1981; Barton and Harvey, 2000). Inserting all the above into eq. (3.1) 
we arrive at the value of A pa 5.0c 2a . Thus the scaling relation for the cortical area 
connectedness in the case when (3 = 0.26 takes the from: 
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Q « 1 - exp (-5.0c 2q \/"' 5 ) , 



(3.3) 



where the exponents a = 0.4, <5 = 0.28, and V g is the gray matter volume of one 
hemisphere expressed in mm 3 . As an example, for the mouse cortex with V g = 56 mm 3 
(Braitenberg and Schiiz, 1998) and since for rodents c ^ 0.15 (Stephan et al, 1981), we 
obtain Q = 0.30. For the human cortex with V g = 3.4 • 10 5 mm 3 and c ~ 0.5 (Stephan 
et al, 1981), we obtain Q = 0.08. Although the brains of these two species differ by 
almost 4 orders of magnitude, their area connectedness differs only by a factor of 3.8. 

Another characteristic of the cerebral cortex related to its connectivity is the so- 
called degree of separation between cortical areas. This quantity is defined as an average 
minimal number of steps (i.e. number of intermediate areas) which are necessary to 
connect two given areas. The detailed calculation of this quantity is presented in the 
Appendix. Here, we only note that the average degree of separation between areas is 
slightly less or around 2 regardless of brain size. 

4 Discussion 

One of the main objectives brains try to accomplish, is to process various forms of 
information. This goal can be achieved by possessing a certain level of computational 
power, which can be controlled by many factors. Some of these factors may involve 
molecular processes, some cellular, and some may involve interactions with the envi- 
ronment. To study all of them, it is an extremely complicated task. However, if one 
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focuses only on a coarse-grained macroscopic outcome of these processes, then a few 
functional/architectonic principles can be identified which brains should meet to effi- 
ciently process information. We suggest that there are three basic principles, which one 
can view as constraints imposed by some hypothetical "perfect" design. We will argue 
below that deviations from this design in real brains are the consequence of compromises 
that brains have to face. 

The first principle is that the number of cortical areas should increase as quickly as 
possible with brain size (Kaas, 1995; Kaas, 2000). The result of Changizi (Changizi, 
2001) on the number of areas vs brain size is consistent with this assertion. That 
trend would allow bigger brains to perform many sophisticated tasks in appropriately 
specialized locations, i.e. locally. This principle implicitly assumes that the number of 
areas is a major contribution to animal's capabilities (Kaas, 1995). The second principle 
is that brains try to maintain a constant connectivity between cortical areas regardless 
of brain size; this is what the results of Sec. 3 may suggest. The third principle, which 
is closely related to the second, is that intra- and inter-hemispheric temporal delay 
should not increase with the brain volume. The last two requirements prevent isolation 
of cortical areas and additionally enhance the efficiency of information transfer between 
them, thus providing the link between local and global information processing (Sporns 
et al, 2000). Although, it has been argued before that delays do increase with brain size 
(Ringo et al, 1994), we treat the third principle as a first order approximation of what 
one may naively expect to be a perfect design. The reason for this is that it is intuitively 
natural to assume that delays would somehow interfere with efficient cross-talk between 
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areas. 

It is interesting to realize that by fulfilling the functional/architectonic requirements 
brains would have to deal with an excessive increase in size of white matter in rela- 
tion to gray matter. This is undesirable, because this would lead indirectly to longer 
cortico-cortical axons, and that would contradict the principle of minimal axon length 
(Mitchison, 1992; Cajal, 1995; Cherniak, 1995; Murre and Sturdy, 1995; Chklovskii and 
Stevens, 1999). Below we analyze how this excessive scaling of white matter vs. gray 
matter arises. 

If we assume that the cortical white matter is composed primarily of cortico-cortical 
fibers then the white matter volume V w is given by proportionality 

V w ~ NL d 2 , (4.1) 

where N is the number of neurons in the cortex and d is the axonal diameter in white 
matter. There has been some debate about whether actually the axonal diameter in 
white matter changes with brain size. The early comparison of fibers for mouse and 
macaque (Jerison, 1991; Schiiz and Preissl, 1996) indicated that majority of axons in 
both species have roughly the same diameter, with the exception of a small fraction 
of thick fibers which were present only in macaque. This led to the impression that 
the fiber diameter can depend only very weakly on brain size, however these data were 
not sufficient to determine the scaling exponent. Recently, the scaling exponent has 
been determined for the group of 6 mammals with volumes spanning almost 3 orders 
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of magnitude using a single methodology (Olivares et al, 2001). It was found that the 
average diameter of axons in corpus callosum scales against the gray matter volume as 
d ~ J/ 0-066 , which indeed is quite weak. 

We want to express the white matter volume in terms of 3 functional parameters: the 
number of areas K, the average connectivity between areas Q, and the conduction delay 
r between areas. For this we have to find how L and d depend on these parameters. 
L is related to K and Q via eq. (2.3). If we expand the exponent in that equation 
and retain only the leading order term, then Q ~ (Lq/K) 2 . Thus, approximately 
L ~ KQ 1 / 2 . The conduction delay r between areas scales as r ~ L /v, where v is 
the velocity of signal transmission in the white matter axons. Since, most of the axons 
in the white matter are myelinated, the velocity v is proportional to the first power of 
axonal diameter d, that is, v ~ d (Hursh, 1939; Rushton, 1951). Thus, the conduction 
delay r ~ L /d. From which, it follows that d ~ L /t ~ KQ 1 / 2 jr. Finally, from the 
Introduction, we have that the number of neurons iV ~ K, ' 9 - Combining all the above, 
we find that the ratio of volumes of white and gray matters depends on the functional 
parameters in the following way: 

V w /V g ~ V 1 ^^- ( 4 -2) 

This equation shows that, if the above three hypothetical functional principles were 
satisfied then the white matter volume would have to grow excessively with the gray 
matter volume due to fast growth of K 3 with brain size. This would lead to an unde- 
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sirable situation when, above certain brain scale, fibers require more space and possibly 
relatively more biochemical resources than the units processing information. A more 
optimal situation would be to require slower increase of K with brain size and simul- 
taneously to allow t to increase, and Q to decrease slowly with brain size, because this 
would reduce a scaling exponent of the right-hand side of eq. (4.2) with V g . That would 
probably slightly decrease information-processing capabilities of the brain but on the 
other hand, it would also decrease its size and biochemical costs. Such strategy seems 
to be taken by brains, since approximately r ~ L /d ~ V®' 20 and Q ~ V~ ' 28 (if we 
take the exponent (3 = 0.26). Thus, in order to manage their biochemical resources and 
size brains probably have to sacrifice some of their computational power by choosing 
compromised scaling exponents. Note that we reached a similar conclusion regarding 
conduction delays as Ringo et al (Ringo et al, 1994), but using a different perspective. 

Another constraint imposed on brains, indirectly related to the above principles, is 
associated with metabolic processes in neurons during information transfer (Laughlin et 
al, 1998). One can expect that for an efficient coding the metabolic energy consumption 
should be low (Levy and Baxter, 1996; Karbowski, 2001; Karbowski, 2002). Below, we 
present a mathematical argument which is consistent with this hypothesis. 

Imaging techniques have been pointing out that glutamatergic excitatory synapses 
are the major users of metabolic energy in the brain (Wong- Riley, 1989; Sibson et al, 
1998; Shulman and Rothman, 1998). Recent estimate of the distribution of energy 
expenditure among different processes confirms that conjecture, especially for bigger 
brains (Attwell and Laughlin, 2001). Thus, the metabolic energy E used by gray 
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matter should be roughly proportional to a number of active excitatory synapses in its 
volume. Since the proportion of excitatory and inhibitory synapses in the cortex seems 
to be brain size independent (Schiiz and Demianenko, 1995), one can conclude that the 
number of excitatory synapses should be proportional to the total number of synapses 
NM in the cortex regardless of brain size. This implies that 

E ~ fNM, (4.3) 

where / is the fraction of active excitatory synapses. Our goal is to find how the fraction 
of active synapses at resting conditions scales with brain size. In order to do this, we 
have first to determine how E scales with the brain volume. From Hofman (1983), we 
find that a basal metabolic rate (an oxygen consumption by the whole cortex at rest) 
E ~ VgVfody 15 ! where Vbody is the body volume. Using a scaling relation between the 
brain and body volumes V g ~ V b o J^ (Hofman, 1983; Allman, 1999), we obtain E ~ V^ 8 . 
On the other hand, the total number of synapses in gray matter, i.e. the factor NM in 
eq. (4.3) scales in first power with V g (Schiiz and Demianenko, 1995; Braitenberg and 
Schiiz, 1998). From this, it follows that the fraction of active excitatory synapses / at 
resting conditions scales with brain size as / ~ V g ~ ' 2 , i- e there is a slow decay in value 
of / as brain size increases. This result can suggest that metabolic energy is indeed the 
resource that brains try not to overuse. 
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5 Appendix 



The degree of separation in a network composed of some units is defined as the minimal 
number of steps which are necessary to connect two arbitrary units. The derivation 
of this quantity for cortical areas is technically similar to the derivation of the average 
degree of separation for neurons, which has been done elsewhere (Karbowski, 2001). 
Here, we follow that approach. 

A given two areas can be connected indirectly through many paths composed of 
chains of other areas in the cortex. The probability S n that these two areas are con- 
nected via at least one of the paths that uses n steps (i.e. the chain of n intermediate 
areas) is given by (compare, Karbowski, 2001) 

S n = l-(l-Q n ) mn , (5.1) 

where m n = (^ K ~ 2 ^jn\ m K 11 " 1 is the number of all possible paths that use exactly n 
steps through which the two areas can be connected indirectly, and Q is the average 
probability of a direct connection between cortical areas given by eq. (2.3) in the main 
text. We can rewrite this formula in a more convenient form as we did in Sec. 2: 

5 n = l-e- p ", (5.2) 

where 
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p n = -K n - 1 hi(l-Q n ). 



(5.3) 



For large n we can expand the logarithm obtaining p n m (KQY^Q. Since KQ is the 
average number of areas a cortical area is connected to, it is substantially greater than 
1. As a consequence p n > 1 as n > 1. We will use this fact later. 

It is also instructive to determine how p n scales with the brain size. From eq. 
(3.2) in the main text, we have to leading order in the Taylor expansion Q ~ V~ s 
(5 = 0.28), and using another scaling of the number of areas K with V g , we obtain 
p n m K n ~ l Q n ~ v g a(n ' 1) ~ 5n = V^ n - a ( n+1 \ Thus p n scales against V g with the exponent 
2(3n — a(n + 1). This implies that p n grows quickly with V g for n > a/ (2(3 — a), and 
decays to zero for n < a/ '(2(3 — a). In the former case, from eq. (5.2) we get 5„ h-> 1 
as V g i— > oo, i.e. there is a connection after n steps in the limit of very big brains. In 
the latter case, this leads to S n i— > as V g i— > oo, i.e. there is no connection at all. 
One can think about the boundary value a/ (2(3 — a) w 1.54 as the average degree of 
separation between cortical areas in the limit of very big brains. As a comparison, the 
average degree of separation between neurons is either 5, for smooth brains, or 10-11, 
for convoluted ones, in this limit (Karbowski, 2001; Karbowski, 2002). 

Because we want to find the minimal number of steps which are needed to connect 
indirectly two given areas, we have to define a second probability P n that areas are 
connected in at least n steps, i.e. n is the minimal number of necessary steps. These 
probabilities are related to the probabilities {S n } in the following way 



23 



P n = (l-S 1 )...(l-S n - 1 )S n 



(5.4) 



for n > 2, and Pi = Si for n — 1. The average minimal number of steps, or equivalently, 
the average degree of separation between cortical areas, defined as < n >= J2n=i n Pn, 
is given by 

K-2 j K-l 

<n>=2-Q+ exp(-^Pi) -(X-l)exp(- Pi)- ( 5 - 5 ) 

j=2 i=l i=l 

This exact but complicated equation can be approximated if we use the fact that p n 
rapidly becomes large as n increases for a given V g . This allows us to neglect all the 
exponents except the leading one, exp[— {p 1 +^2)]- Thus, approximately 

< n >^ 2-Q + (1 -Q)e~ KQ \ (5.6) 

and it is apparent that the average degree of separation between areas is close to 2 with 
only minor species specific corrections and finite size effects. For example, for rat with 
K = 26, Q = 0.25 we find < n >= 1.89, for cat with K = 65, Q = 0.27 we have 
< n >= 1.73, and for macaque monkey with K = 73, Q = 0.15 we obtain < n >= 2.01. 
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Figure Captions 
Fig. 1 

Dependence of the average size of cortical patches in VI on the logarithm of the gray- 
matter volume V g for several mammals. Data for the patches widths are: for rat £ = 370 
(Km (Rumberger et al 2001), for tree shrew £ = 230 /im (Rockland et al 1982), for cat 
£ = 300 pim (Luhmann et al 1986), for macaque £ = 230 fxm (Amir et al 1993), for 
human £ = 400 /im (Burkhalter and Bernardo 1989). The gray matter volumes were 
taken form Stephan et al, 1981. 

Fig. 2 

Log-log dependence of the average axon length L Q in white matter on the gray matter 
volume V g in one hemisphere for few mammals: mouse, rat, cat, and macaque. We 
find the allometric law L = 1.10 Vg . The scaling exponent (3 = 0.26 ± 0.02. Data 
for the axon length come from computation performed in Sec. 2. For macaque, it was 
taken L = 14.5 mm as an average value between 12.0 and 17.0 mm. The gray matter 
volumes were taken form Stephan et al (1981) and rescaled by 0.5. 



32 





1 1 


i i i i 




rat 




human 



cat 




tree shrew 




macaque 












i i 


i i i i 



log V g (mm" 3 ) 



Figure 1, Karbowski 



